Solitons of Bose-Fermi mixtures in a strongly elongated trap. 
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It is shown that a Bose-Fermi mixture of a degenerate gas of spin-polarized fermions, whose 
number significantly exceeds the number of bosons, embedded in a strongly anisotropic trap, is 
described by the one-dimensional coupled nonlinear Schrodinger equation for the boson component 
and the wave equation with external source for the fermion component. Depending on the type 
of boson-fermion interaction, the system may display modulational instability and the existence 
of solitons in the fermion and boson components respectively. Such solitons represent either a 
local decrease (increase) of the density of both the components or a decrease of the density in one 
component and an increase of the density in the other component. It is shown that the type of the 
effective interactions can be easily managed by varying the trap geometry or by means of Feshbach 
resonance. 
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I. INTRODUCTION 



Stable coherent structures represent one of the main 
subjects of interest in the mean-field theory of quantum 
gases at zerotemperature [1]. In quasi-one-dimensional 
condensates of bosonic atoms, such structures are con- 
ventionallyassociated with solitons: dark solitons in the 
case of repulsive interactions among atoms, and bright 
solitons in the caseof attractive interactions. Successful 
experimental generation of dark and bright solitons has 
been reported inBose-Einstein condensates in Refs. [2] 
and Refs [3], respectively. Considerable progress has al- 
sooccurred in the experimental creation of more sophis- 
ticated condensates, in particular, boson-fermion(BF) 
mixtures (see e.g. [4, 5]). The first theoretical studies of 
such systems have revealed anumber of interesting prop- 
erties. One such property is an effective attraction that 
develops among bosonsdue to the boson-fermion inter- 
action [6] which affects the background distribution and 
the stability properties of thecondensate [7]. Very re- 
cently, numerical work [8] has shown that localized dis- 
tributions ofcondensed atoms, associated with solitonic 
trains, can form in a strongly anisotropic trap even when 
the boson-boson interaction is repulsive. In Ref. [9] it 
has been suggested that the respective localized states, 
can be interpreted as a matter wave realization of quan- 
tum dots. A boson-fermion mixture can exhibit quite dif- 
ferent types of behavior according to the relative ratio 
of the two components. The model explored in [8] and 
[9] addresses a mixture of a relatively small number of 



fermions embedded in a much larger bosonic component. 
In this paper, we consider the opposite situation of a rel- 
atively smaller number of bosons embedded in a large 
fermionic component. Such a large fermion component 
can become degenerate, and only a small portion of the 
fermions in the vicinity of the Fermi surface participates 
in the kinetic processes. If the fermions were also spin- 
polarized, and thus noninteracting, they would represent 
a linear system. It turns out however [6], that the pres- 
ence of bosons, even in a relatively small number, can 
introduce nonlinearity into the fermion system, thereby 
significantly changing the properties of the system. In 
the present paper we obtain dynamical equations gov- 
erning the effectively one-dimensional (1-D) BP mixture, 
show, in particular, the possibility of the appearance of 
two-component (BF) solitons in a strongly anisotropic 
trap, and describe some of the simplest properties of the 
solitonic solutions. It is to be noted in this context that 
two natural limits exists for bosons embedded in a Fermi 
cloud, depending on the relative densities of the fermions 
and bosons. One limit is when the boson number domi- 
nates the fermion number and is the subject of the anal- 
ysis in Refs. [8, 9]. The other limit is when the fermion 
number dominates the boson number. The latter case, as 
well as a discussion of the differences in the two cases, is 
considered in the present paper. The paper is organized 
as follows. In Section II, we describe the "Bose-Fermi" 
system. Section III explains the various length scales of 
the system. Section IV shows the reduction of coupled 
3D equations to a single 1-D equation that governs the 
system. In Section V we summarize our main results and 
discuss their significance. 
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II. STATEMENT OF THE PROBLEM 



malization conditions are 



A mixture of bosons and spin-polarized fermions, with 

a dominant fcrmionic component, can be described by 
coupled mean-field equations for the boson order pa- 
rameter ^(r,f)and the fermion densityn(r, = no(r) -|- 

6n{T,t), where no(r) = ^ [^{Ef - ^^(r))]'/', which 
were derived in [6]: 

+ TTItA* - Vb^ - 511*1* - 52"* = (1) 



dt 2M 



noV 



1/3 



m 



(2) 



Here M and m are the respective masses of bosons and 
fermions, the two-body interactions are characterized by 
the coefficients.gi = , and 52 = ^2^^ with 

scattering cross-sections Ubb and a^f, and /x = 
(see e.g. [7]), no(r) is the fermion density described 
by the Thomas-Fermi (TF) approximation for degen- 
erate Fermi gas with Fermi energy Ep, and 6n(r,t) is 
the departure of the fermion density from no(r), con- 
sidered to be small enough: |(5n(r,t)| ^ 'T.o(r)- We 
consider the trap potentials V for both components 
to be parabolic and strongly elongated along the x- 
directioniV^B = M^{A'^x^ + )and Vp = ^{X^x^ + 
r^) .Here A and A ^ 1 in the boson and fermion 
potentials Vb and Vp respectively, represent the high 
anisotropy, i.e., the respective aspect ratios of the traps: 
n and Lo are the linear oscillator frequencies of bosons 
and fermions in the transversedirection, and r± = {y, z). 
Introducing the dimensionless variables T = ifit and 
R = ^, where a = sJhjM^ is the transverse linear oscil- 
lator length of the bosons, and re-scaling through 



V(R,r) 

p(R,T) = 



= 2^7r|a66|a*(r,t); 
47ra^|a6/|n(r,i); 



(3) 



allow the TF distribution for thefermion density no(r) to 
be written as 



where Kp 



po = K.F {Sf - (A^X^ + R\)\ 
3 



3/2 



(4) 



2 (^)^ andf^ = 2M£2|£ is the di- 
mensionless Fermi energy. Then Eqs.(l), (2) acquire the 
form 



+ ARV-(A2x2+iii)V>-2ai|^|V 

- (72 pip = 



sVi 



Vr 



PoVr 



\ 1/3 



Pi + X0-2|V'P 



(5) 
(6) 



where Pi(R,T) = /j(R,T) - po(R), 0-1,2 =sign(5i,2), 



a 



321 (M.Y 

2 V m/ 



2/3 



andx 



abf_ 
abb 



.The nor- 



j IVl'rfR = 47r^7V6 and J pdR = 4^^^/ , (7) 
Nb and Nf being the total boson and fermion numbers. 



III. SCALING 

In the present paper we consider the case when the 
transverse dimension of the Fermi component, Rp, sig- 
nificantly exceeds the transverse dimension of the Bose 
component, Rb- The orders of the respective scales can 
be identified as Rp — \fEp [which follows from (4)] and 
Rb = 5, which in dimensionless units corresponds to 
the transverse oscillator length of the low density Bose 
component. Thus we have, £p » 1. Multiple-scale ex- 
pansions [10], involving a small parameter identified as 
the ratio between the transverse oscillator length and 
the healing length, provide a self-consistent and control- 
lable approach in the case of a one-component BEC, 
and reduce the 3D Gross-Pitaevskii equation to an ef- 
fective ID model [11]. The use of ID models, in turn, 
reveal interesting features of the low-dimensional dynam- 
ics, allows one to obtain approximate analytical solu- 
tions, and results in dramatical reduction of computa- 
tional resources in studies of the evolution processes. 
Since BF mixtures are two-component systems charac- 
terized by a variety of scales, a natural question arises: 
Is such a reduction possible in a Bose-Fermi mixture! 
This is the question we address in our present study 
through a generalization of the approach of ref. [11], when 
the fermion number significantly exceeds the number of 
bosons. To this end we first identify the small parame- 
ter e = I <C 1, where ^ = (87rnh|a(,b|)^^^^ is the boson 
healing length (rib is the maximal boson density). We 
restrict oiu analysis to fermion numbers large enough to 
satisfy Rj^/R^ ^ £p = ep /e'^ , where ef ^ 1 is a param- 
eter of the problem. Then we can describe the condition 
£p 3> 1, using the same single small parameter e. This 
scaling allows us (by analogy with Ref. [6]) to concen- 
trate on the spatial domain, limited in the transverse 
direction by some radius R, which satisfies the condition 
Rb "C -R <C Rp. Physically, such scaling is available due 
to Pauli's exclusion principle which forces the fermions 
to occupy a much larger region relative to the bosons. 



Specifically, we impose R 



-1/2 



The elongation of the 



fermion trap is expressed through A < e^^^. The spa- 
tial domain of the fermions in the longitudinal direction 
is defined by the constraint, e^X^X^ <C 1, which gives 
^ ^ while for the case of bosons when A ~ A ~ e^/^, 
we have that X < e'^/^ Thus we consider the situation 
in which the boson component is embedded inside the 
fermion system. Subject to the above conditions the TF 
distribution (4) can be expanded in the Taylor series 



Po 



: tCpSl/'^ 



l-^^X^X^ + Rl)+0{e^) 
2ep 



(8) 
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We look for a solution in the form of the seriesi/'(R, T) = 
eVi + £^-02 + • • • , andp(R, T) = pq + e^pi + • • • .To justify 
the above expansion, we introduce a characteristic width 
I of the bosonic excitation in the axial direction so that 
the boson density may be connected to their total number 
by means of the estimate: Nh ~ wla^n. This allows 
us to rewrite the small parameter as e ^ ^/8N\abb\/i. 
In the case of a localized distribution, £ ~ ^, while for 
the estimates involving plane waves and dark solitons 
£ ^ L where L = a/\fh. is the longitudinaldimension 
of the condensate. Then, for excitations with ^ ^ ^, 
from (7) we readily obtain the estimate li/^ip ^ l.On 
the other hand, the scaling we have introduced implies 
IV'P ^ Ip^PoI which corroborates the fact that the boson 
and fermion numbers involved in the dynamics are of the 
same order (see also [6]). Indeed in momentum space, 
the fermions in question are located in the vicinity of the 
Fermi surface; the actual number of fermions, p, is much 
bigger than its time dependent part pi. In other words 
\p — Po\ < Po is verified. 



IV. ONE-DIMENSIONAL EQUATIONS 

The natural occurrence of the small parameter e in the 

problem suggests a perturbation treatment in powers of 
e. A regular analysis, however, gives rise to secular terms 
and, consequently, to an inaccurate description. Hence 
an application of multiple scale analysis becomes neces- 
sary. Following standard procedures (see e.g. [11]), we 
introduce scaled variables tj = e^T and Xj = e^X, re- 
garded as independent variables. Substituting the above 
expansions for time and space into eq. (5) and gathering 
the terms of the same order in e wc rewrite the equation 
for a boson macroscopic wave function in (1), in the form 
of a set of equations: 
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Here the respective eigenfunctions of the operators are 
defined as, C±Cim = £imCim{'R±) and£a,(?!)„(a;o) = 
Vn^Pnixo) where rj^ = (2n -|- l)z/^, = 2{l + m. + l)iy±, 
and the scaling chosen guarantees that u± > O.The ba- 
sic idea of our analysis is to look for a solution to (5), 
which describes the evolution of the background state of 
the mixture, in the form: 

^1 = v-'/*A{xuh)MxoKo{B.i^)e-"'°'° • (11) 

The envelope function in eq. (11) is a function of slow 
variables, {xi,X2,-" ,ii,i2,---)) from which only the 
most rapid are indicated explicitly, and is independent of 

xo, to andi?_L.Ineq. (11) (f>o{xo) = (^)^^^exp {-l^xXo) 

1/2 

and Co{R±) = (^) exp (— ij/^i?^) arc the normal- 
ized eigenfunctions of £x and £±, i.e., the ground state 
wave functions of the ID and 2-D harmonic oscillators, 
respectively. The factor Vx is introduced to pro- 
vide the requirement |Ap ~ 1 so that the condition 
\tpi\^ ^ 1, established above is satisfied. By direct 
substitution one verifies that ansatz (11) satisfies the 
equation obtained in the first order ofe (j = 1) when 

'^0 = <^2^^f£'p'^ — {i^x + 2z^_l). The primary goal is 
the determination of the time evolution of the envelope 
function, A{xi,ti) as a solvability condition [10] for the 
multiple scales analysis. To this end we pass to the sec- 
ond order {j = 2) equation, and find that 



dtpi 



^2 = —i-^T' 



dti dxodxi 



(12) 



Using an ansatz for ip2, such that it is given as the 
expansion over the complete set of the eigenfunctions 
(f>n{x) , Cm{R±) and subjecting it to the intial condition, 
that there are no excited modes Cm{R±) with m > 0, '^2 
can be written as, 



(9) 



V'2 



-1/4 



Bn<PnCo exp{-iu]oto), (13) 



where j = 1, 2, Fi = 0, 



dti dxodxi 

F = i^'^^ 2 ^^^^ ^'^^^ i^^^ 2 ^^^''^ 
9*2 dxodx2 dx\ dti dxodxi 

M 

-|-2(7l|V'l|^'!/'l + 0'2 Pli'l 

M 

and j = 1, 2, ... coincides with the order of e at which the 
equation is obtained. The operator £ in eq. (9) is given 

by 

£ = £1. + £x ^ 

M 

£^=-/^^ + vlRl, ul = l-a2^)CF£'J^ (10) 



~ a^2 + ^x^O 



Q^2 ' "x-O, "x-^^ (72— /VfC^ 



Substituting the above form of ip2 into the second order 
equation (9), we find, 

n ^ ^ 

We can see from the above condition that the amplitude 

Bn becomes infinite, when n = 0. Such terms arc referred 
to as "secular terms" [10]. The method of multiple scales 
allows us to identify and eliminate such secular terms by 
requiring F2 to be orthogonal to ipi. From this condi- 
tion, wc can obtain the solvability condition for A to be 
OA/ dti — and thus A does not depend on ti. We write 
A — A(f2-,xi) and when n ^ 0, we find the solution to 
the second order. 



= ,.-1/4 



i>2 = v.. 



OA 



- n.O 



9XI VO - Vn 



^{xo)Co{-R^)e-'^°'". (14) 
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where r„^o = -2 J^^4>n{x)-£:(j)o{x)dx.The equation of 
the third order, however, involves the dependence on the 
fermion density, which means that it must be solvedto- 
gether with Eq. (6), which in the new, scaled variables, 
after substituting for tpi from (11), acquires the form 



dtl 



■v%Vl^pi+x{Ro,to)\A\' 



where Rq = {xo,y, z) and Vp = aJC^^Sp 



2 ^2 



(15) 



(16) 



arid 

3 2m '^-nCl 



In the physical variables, v|. = 

therefore in what follows, we refer to Vp as the Fermi 
velocity. Now Eq. (15) acquires a transparent physical 
meaning. Fermions as quasi-particles are excited on the 
Fermi surface, and thus propagate with the Fermi ve- 
locity along the condensate. Their dynamics is described 
by the wave equation with a source term, which describes 
the effect of the boson distribution on the fermion den- 
sity.Because A{xi,t2) is independent of R±, xq and to, 
it can be treated as a constant in Eq. (15). Since (f>Q(xo) 
changeson the scale of the longitudinal dimension of the 
boson component, £ ^ ^, while (^o{R±) changes on the 
scale of transversedimension a, the effect defined by the 
curvature of (q{R±) dominates in x'and we can neglect 
the second term in Eq. (16). We are interested in cou- 
pled fermion-boson excitations (in particular, solitons) 
when pi is represented by the particular solution of the 
in homogeneous equation (15) as shown below: 



Comparing pi to ipi from eq. (11), given by 



TT 



(17) 



(18) 



we find that pi is also governed by the same envelope 
function A{xi,t2) as the bosons. Substituting this ex- 
pression for pi in F3 we can write the third order expan- 
sion {j = 3) for the bosons as 



OT2 



^ dtp2 _ ^ 

dti dxodxi 



dxodx2 dxf 
+ 2ai|Vi|'V'i 



(7ra)3/2 p 



lApe-'^-R^Vi- (19) 



Here we have used the fact that a| = 1. Requiring the 
absence of secular terms in the third order equation of the 
multiple scale expansion, which is the same as requiring 
the orthogonality of F3 with ^2 and ipi , and employing 
explict expressions for tpi and ip2, given by (18) and (14), 
respectively, we find that A{xi,t2), satisfies. 



.dA 



d'^A 



0, 



(20) 



which is the solvability condition we seek. In eq. (20), 

^ = 1 + Er=i £i = 1/2; Xeff = -X66 + Xbf, is the 
effective nonlinearity characterizing the combined effect 
of: direct boson-boson interactions, Xbb = -j^^^v^, and 
the interactions among bosons induced by boson-fermion 



-v^. The solution to 



scattering processes, XhS - 2^3/2^3/2^ ^ 
equation (20) is known exactly. This is the main result of 
our analysis and shows that bright and dark soliton so- 
lutions are possible for the envelope function. According 
to the method of multiple scales, t2 = e^T and x\ ^ eX , 
thus we can re- write the solution, A{xi,t2) of the eq. (20) 
in terms of our original scaled variables, X and T, as 
A{eX,e'^T). Eq. (20) along with eqs. (18) and (17) de- 
scribe the evolution of the original system as defined by 
(5) and (6) for times up to 1/e^ (e ^ l).When the boson- 
fermion interaction is turned off (xbf = 0), we expect 
the boson and fermion components to be independent of 
each other, i.e., eqs. (5) and (6) are decoupled. In such 
a case, Xeff = ^Xbt, and from eqs. (17), (18) and (20), 
we see that the fermion density is given by the Thomas- 
Fermi Approximation, eq. (4) and the boson system is 
described by NLS equation as expected. When Xbf 7^ 0, 
we see that boson-fermion interactions always lead to at- 
tractive (effective) interactions between the bosons. This 
qualitatively agrees with earlier results [6]. Our analysis 
below shows that the quantitative dependence of the be- 
havior of the mixture on the system parameters is very 
different for the case when the number of fermions is 
much less than the number of bosons compared with the 
case when the number of fermions is much larger than 
the number of bosons. 



V. DISCUSSION AND ESTIMATES FOR REAL 
CONDENSATES 

Since our analysis has resulted in the NLS equation 

(20) governing the dynamicsof small amplitude excita- 
tion in BE mixtures, we know at once that several well- 
known consequences of the NLS equation will follow, in- 
cluding modulational instability of plane waves, and exis- 
tence of bright and dark solitons, depending on the sign 
of Xe//- When Dxeff > 0, the envelope function, A 
is given by the bright soliton solution ("sech(eX)") and 
when Dxeff < 0, it is given by the dark soliton solution 
("tanh(eX)"). Since D > 0, the bright or dark soliton 
solution is solely dictated by the sign of Xef f ■ In the case 
where boson-boson and boson-fermion interactions have 
different signs, the observable phenomena depends on the 
relative values of Xbb and Xbf- One can easily find the 
critical value for the boson-fermion interaction at which 
both interactions are balanced, i.e. when Xe// = 0: 



Ml 



MEf 



1/2 



y/a abb 



(21) 



This value that we obtain differs from that reported in [8] . 
This is expected given that the underlying physical con- 
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ditions are different in the two cases. Note, in particular, 
the weak dependence of a^J in eq. (21) on the fermion 

number, |agj| ~ N^^^^"^ in contrast to the dependence 

\^bf\ -^F^^^ reported in [8]. This is due to the fact 
that in our case we consider a small number of bosons 
interacting with a large system of degenerate Fermi liq- 
uid and changing the number of fermions does not have 
much impact while in [8] the small number of fermions are 
trapped by a large system of bosons. In such a case, the 
effect of interaction with each fermion becomes impor- 
tant and hence a large dependence on the fermion num- 
ber. In our result, we also find that |agy| ~ ^/ahb while 

in [8] |ag^| ^ a^^^.When \abf\ equals the critical value, 
the system becomes an ideal gas, described by the lin- 
ear Schrodinger equation. If \abf\ is larger than the above 
critical value, Xbf > Xbt, interaction between bosons and 
fermions dominates, leading to a negative effective attrac- 
tive interaction between the bosons. This is the case 
where modulational instability and bright solitons can be 
observed. Eq. (21) shows that value of |agj| can be ma- 
nipulated by changing the radial trap size, o.Note that 
using the parameters of the Rb^^-K^" mixture experi- 
ment [4]: n = 27rx 197 Hz, a^b = -1-5.25 nm,and the Fermi 
temperature 0.25 jjK, allow us to compute |agy | « 32 nm. 
On the other hand, using the data of Li^-Li*' mixture, 
Ref. [5]: Q = 16.14 kHz, Ubb = -1.35 nm, the Fermi 
temperature 4 jjK, gives us |aj;j| « 2.32 nm. An alter- 
native way of making the boson-fermion interactions to 
be dominant is to increase abf (or decrease |af,(,|)[8].An 
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FIG. 1: Two-component soliton formation along the elon- 
gated trap direction (x-dircction) for the Rb^^-K'"' mixture 
example. The solid line shows the density of bosons, |^'(r,t)|'^, 
corresponding to the bright soliton solution for the envelope 
function. The dash and dash-dot lines shows the fermion den- 
sity, 5n{r,t), for repulsive and attractive interactions between 
bosons and fermions, respectively. Taking into account the 
nonzero background, the excitation of fermions in the former 
case can be associated with a conventional dark soliton, while 
in the later case, it is a bright-against-a-background soliton. 

interesting peculiarity of the bright soliton solution of 
eq. (20), follows from (17). If the fermion-boson inter- 



action is repulsive(cr2 > 0), the fermion density, has a 
relative negative sign compared to the boson density and 
is lower in the region of the soliton than outside. In 
thiscontext, the terminology "bright soliton" does not 
reflect the situation precisely, because in the fermionic 
component oneobserves darkening. To illustrate this, we 
plot the densities of bosons (solid line) and fermions for 
the for (72 < (dash-dot line) and fT2 > (dashed line) 
along the elongated trap direction (x-direction) in Fig. 
1 for a Rb^^-K***^ mixture [4], when the envelope solu- 
tion of eq. (20) is a bright soliton (Xeff > 0). While 
for attractive boson-fermion interactions one haseither 
bright-bright or dark-dark solitons, in the case of repul- 
sive boson-fermion interactions, one can excite bright- 
dark anddark-bright solitons. Here, the first and sec- 
ond term refer to bosons and fermions respectively (see 
also Fig. 1 for the explanation of the terminology for the 
fermions). For the calculations presented in the present 
paper, the choice of an appropriate scaling in finding the 
small parameter e, was of primary importance. In or- 
der to check its feasibility for experimentally available 
atomic gases, we take into account the possibility of a 
change of the scattering length by meansof Feshbach res- 
onance and consider an example using Rb^^-K^" mix- 
ture [4], with a radial size ofa « 2 /im, and longitudinal 
extension L « 200 /xm with the Fermi temperature of or- 
der of 0.25 /xK. Takingthe number of rubidium atoms to 
be 500, the scattering length a^bu = 0.1 nm, and the width 
of the soliton £ ^ 20 /im, we obtain the mean healing 
length C ~ 14 ^m, which gives e = 0.14,5^ = 3910 > 1 
and ep ~ 76.6. We also notice that for the case of the 
Li^-Li^ mixture [5] with the Fermi temperature 4 fxK we 
obtain £p = 500 » 1, which also matches with the con- 
ditions of the application of the theory. Fig 1 shows the 
formation of the two-component soliton for such a Rb^^- 
K^*^ mixture with |at/| = 20 nm, when the boson-fermion 
interaction is attractive (fT2 < 0) andrepulsive (a2 > 0). 



VI. CONCLUSION 

In this paper we have considered a boson-fermion mix- 
ture, where the number of fermions is much larger than 
the number of bosons and is confined by a strongly 
anisotropic trap. It has been shown that both the bo- 
son and fermion systems are described by an effective ID 
nonlinear Schrodinger equation (20) and thus can display 
modulational instability and the existence of bright and 
dark solitons. Such solitons involve most of the bosons of 
the system and a relatively small portion of the fermions 
in the vicinity of the Fermi surface, where the soliton 
propagation is along the axial direction of the conden- 
sate. In a system without coupling between the bosons 
and the fermions, the fermions would tend to propagate 
with their Fermi velocity. However, the coupling to the 
bosons that is present in our system forces them to form 
a slow (and even static) solitary excitation following the 
behavior of the bosons. A remarkable feature of the result 
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obtained is that the soliton dynamics is governed by a sin- 
gle nonlinear Schrodinger (NLS) equation, even though, 
physically, it describes two separate systems. The soli- 
ton solution of the NLS equation can describe either an 
increase or a decrease in the densities of the two compo- 
nents. In the case of the attractive interaction between 
the components, the solitons may describe simultaneous 
increase or decrease of the densities of both the compo- 
nents of the mixture. In the case of the repulsive inter- 
actions between the components on the other hand, they 
can describe an increase of the density in one component 
and a decrease of the density in the other component. 
This is quite an unusual situation: there is only one soli- 
ton solution to the NLS equation for a given effective 
nonlinearity, however it describes opposite effects in the 
two components. Another feature to be mentioned is that 
due to the finite extension of the condensate in the longi- 
tudinal direction, an internal mode is excited. This exci- 
tation results in the effective mass of the bosons, making 
up the soliton, being twice as large as the real mass, 
which does not happen in the case of infinite extent of 



the trap. We have explicitly shown that the type of the 
effective interactions can be experimentally manipulated 
either by just affecting the trap geometry or by means 
of Feshbach resonance (the latter option being also effi- 
cient in the opposite case of relatively small number of 
fermions [8]). We hope that observational attempts to 
verify these physical aspects of the results we have ob- 
tained will be made in future experiments. 
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